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Abstract 

Using bases of fourth rank tensorial bases of [[V 2 ] 2 ] symmetry elaborated by Walpole we obtained 
expressions for inverse of Young's modulus E(n), inverse of shear modulus G(m, n) and Poisson's 
ratio i/(m, n), which depend on components of the stiffness tensor S, on direction cosines of vectors 
n of uniaxial load and the vector m of lateral strain with crystalline symmetry axes. Crystalline 
media of high and medium symmetries are considered. Such representation yields decomposition 
of the above elastic characteristics to isotropic and anisotropic parts. Expressions for Poisson's 
coefficient are well suited for studying the property of auxeticity. 
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I. INTRODUCTION 



Recent interest in single crystal materials has made the understanding their elastic char- 
acteristics, i.e. the compliances Cy and stiffness coefficients [lj], the bulk and Young's 
moduli E, Poisson's ratio v and the shear modulus G, increasingly important j^-j^- The 
anisotropy of Poisson's ratio was studied in several papers In particular there is 

fundamental and practical interest in materials with negative Poisson's ratio, called aux- 
etics Q. Love presents a si„ g le e xaI „p, of Cubi o ^ crysM p y nte as h avi„ g a 
Poisson's ratio of -0.14. Lakes described the synthesis of an actual auxetic material and 
proposed a simple mechanism underlying the negative Poisson's ratio Q|. Alderson and 
Evans fabricated microporous polyethylene having a negative Poisson's ratio In spite of 
general belief it was shown that many crystalline solids, among them metals, have a negative 
Poisson's ratio Ql-Q]. Data were compiled from the literature by Baughman et al. and 
were analyzed to show negative Poisson's ratio to occur for stress in an oblique direction 



upon single crystals of cu 
In our previous paper 



)ic metals 



2l|, we considered Poisson's ratio of all cubic materials stretched 



in [001] direction and measured for [100] lateral direction. We concluded that auxeticity is 
a rather common phenomenon. 

Studying auxeticity of cubic elastic media we noticed that the size of the auxetic part 
of the stability region depends on the direction n of stretch and the direction m of lateral 
strain. To verify this observation we derived general expressions for E(n), u(m, n), and 
G(m, n) for all stable cubic materials. 

Our derivation gives expressions which do not depend on the choice of Cartesian coor- 
dinate system. To obtain such formulas, we used fourth-rank tensorial bases introduced by 
Walpole 12]. This approach is simpler than the method developed by Turley and Sines 
Q, Q] and applied by Gunton and Saunders [3, 0, 0], Li j^| and Jain and Verma To 
obtain Poisson's ratio and Young's and shear moduli for an arbitrary direction of applied 
stress, these authors used the method of rotations by Eulerian angles. 
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II. POISSON'S RATIO, YOUNG'S MODULUS AND SHEAR MODULUS 



Elastic properties of an anisotropic body are characterized by the elastic constants Ciju 
or by the elastic compliance coefficients Sijki relative to an orthogonal coordinate system. 
Consider any two specified orthogonal unit vectors n and m and three related characteristics 



of elastic properties of the considered body |23j, |24(, namely the Poisson ratio z/(m, n 



i/(m, n) = = , (1) 



Young's modulus E(n) 



E(n) 



and the shear modulus G(m, n) 



riirijSijkirikni, (2) 



miTijSijkimkni. (3) 



4G(m, n) 

Repetition of a suffix in a product of tensors or in a single tensor implies the usual summation 
with the respect to that suffix over the values 1,2,3. 

Introduce three n ® n and m <g> m and m <S> n with elements riiUj, rniTrij and rriirij 
(i,j = 1,2,3). Using these dyads the definitions ( fTj) -([3*f can be written in a brief form (cf. 
Rychlewski \2^). 

£(n) -1 = (n®n) • S- (n®n), (4) 

[4G(m, n)]- 1 = [4G(n, m)]- 1 = (m ® n) • S • (m ® n), (5) 
v(m,n) v(n,m) 

■■ — (m (gi m) • S • (n ® n). (6) 



E(n) E(m) 

In an arbitrary chosen coordinate system, both £'~ 1 (n) and [4G(m, n)] _1 are related to 
quadratic forms of a positive definite matrix, hence £'~ 1 (n) > 0, [4G(m, n)] _1 > 0. 
Tensors C and S with elements C^-jy and S^jy, respectively, obey the relation 



or for components 



This means that S = C and 



C • S = S • C = I 4 , (7) 



CijrsS-rskl — Sij rs C rs kl — ($ik$jl ^il^jk) ■ (8) 
-1 
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III. POISSON'S RATIO, YOUNG'S MODULUS AND SHEAR MODULUS FOR 
ISOTROPIC ELASTIC MEDIA 



Introduce two fourth rank tensors - J with components 

Jijki = g (h)ij (h) kl , (9) 

and K = (I 4 — J) (cf. Walpole [ssj])- The second rank tensor I 2 has components (h)^ = 
The introduced tensors obey several important relations 

J 2 = J, K 2 = K, 

J + K = I 4 , J K = K J = 0. (10) 

The tensors J and K have complete Voigt symmetry. Using relations (|7j)- (jl(J|) . one obtains 

C = cjj + ckK, 

S = sj3 + s K K, (11) 

where Q 

cj = ^ = C\\ + 2C12, ck = — (Cijij — -Ciijj \ = {Cxi ~ C12) , sj = Cj 1 , sk = . (12) 

This means that S = cJ 1 J + K. An isotropic medium is mechanically stable if cj > 0, 
c K > {sj >0,s K > 0). 

The tensor S can be rewritten in a useful form 

S = {sj - s K ) J + s K U. (13) 

From the definition (@J) and the form ()13J) of the tensor S, one obtains 

E-\n) = E- 1 = {sj - s K ) /3 + s K = sj + 2s K /3 > 0. (14) 

Hence, using the tensorial basis we obtained the familiar expression for £' _1 (n) 2^|. 

Calculating products (m ® m) • S ■ (n <g> n and (m <g> n) • S • (m <g) n), we obtain 

Km,n) = (sj - sk) _ _Sx2 = C 12 
E 3 3 cjCk 

1 ^ = ^>0. 



4G^,r^_ | 2 2c K 

These formulas can be found in books |2c 

In agreement with Ting and Barnett [26|, an isotropic medium is completely auxetic (i.e. 
f(m, n) < for all pairs (m, n)), if £12 > (or Cxi < 0) and nonauxetic (i.e. i/(m, n) > 
for all pairs (m, n)), if < (or Ci 2 > 0). 



IV. POISSON'S RATIO, YOUNG'S MODULUS AND SHEAR MODULUS FOR 
CUBIC ELASTIC MEDIA 



As the three crystallographic directions of the cubic system are mutually perpendicular, 
the unit vectors a, b and c satisfy the conditions 

aa = bb = cc= 1, ab = ac = bc = 0, a®a + b(g>b + c(g)c = l2. (15) 

Introduce the symmetric fourth rank tensor £ with components 

= aidjdkai + bibjb k bi + CiCjC k ci. (16) 

With the help of S, Walpole constructed two different fourth rank tensors L and M 

L = I 4 — £, M = £ — J. (17) 

These tensors obey the following relations 

L + M = K, J + L + M = I 4 , 
L 2 = L, M 2 = M, 
JL = LJ = 0, JM = MJ = 0, LM = ML = 0. 

For all cubic symmetry classes (23, m3, 432, 43m and m3m) the tensors C, S can be 
expanded in the basis of tensors J, L, M 

C = cjj + c L L + c M M, 

S = sjJ + s L L + s M M, (18) 



where 22 



cj — C\\ + 2Ci2, cl — -LijkiCijki — 2C44, cm — -MijkiCijki — C\\ — C12, 

cj, cl and cm are the eigenvalues of the tensor C. Similarly, sj, sl and sm are eigenvalues 
of S. The mechanical stability is guaranted when all these eigenvalues are positive. 
The tensor S (J 18)) can be written in a useful form 

S = (sj - s M ) J + s L I 4 + (s M - s L ) S. (19) 



The coefficients cj, cl, and cm can be written in terms of Six, Sx2, and 544, namely 



sj = Sxx + 25*12, s L = 544/2, s M = Six - Sx2- (20) 

Using these relations and second of Eqs. (|19|) we express the matrix elements SV,- by CV,- 
(i,3 = 1,2,4) 

(Cxi + 2C12) (C11 — C12) (Cxi + 2C12) (Cn — C12) C44 

These expressions are in agreement with familiar results 25]. 

We calculate Young's modulus in direction n with the help of formulae (JBJ) and l|19|) . As 
a result, we obtain the familiar formulae 

E-\n) = [( SJ - s M ) /3 + s L ] + (s M - s L ) T(n) = £ u - 2 (S n - 5 12 - S u /2) 

/ 9 2 2 2 2 2 \ 

x (™ a % + n fe n c + n c n a ) , 

where T (n) = J2i= a n t- 

Similarly, for z/(m, n) and G(m, n), we get 

v{ m n't 
£(n) 

[4G(m, n)]- 1 = s L /2 + (s M - s L ) P(m, n). 
The function P (m, n) was introduced by Ting and Barnett \i\ 

c 

P (m, n) = P (n, m) = £ (m^) 2 , (23) 

i=a 

where, for example, m a = (ma). Note that the first term of expression defining 
[— z/(m, n)/E(n)] is isotropic, whereas the second one is anisotropic. The same structure 
has the shear modulus. 

The ratios v/E and [4G] _1 are not independent 

viva, n) 1 1 . .1 . . 

- " 77v V = q Sj ~ SM ~ o SL - 24 

E(n) 4G(m, n) 3 2 

If C44 = (Cxx — C12) /2, one deals with an isotropic medium. In this case, sl = sm = $k, 
therefore the last terms of i? _1 (n), u(m, n)/E(n) and [4(7 (m, n)] _1 vanish, and one gets the 
results obtained in Sect. IIHI 
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Using Eq. (|20|) one can write z/(m, m) (|22p in the form obtained by Ting and Barnett 

H 

v(m, n) = [1 - 2P(m, n)] S 12 + QP(m, n)/2, (25) 
where Q=[2 (Sn + S12) — 544]. The functions T(n), P(m, n) obey the inequalities 

1/3 < T(n) < 1, 
< P(m,n) < 1/2. 

The function p(n) achieves the minimal value for n = (111) and maximal value for n = (001). 
The function P(n, m) attains the minimal value for n (or m) equal to (001) and for 
arbitrary unit vector m, perpendicular to this n (or for arbitrary unit vector n perpendicular 
to m = (001)). The maximal value 1/2 of P(n, m) is reached for n equal to (—110) or (110) 
and m = (110) or m = (—110). 

The ratio [— z/(m, n)/P(n)] obeys the inequality 

S 12 < -i/(m, n)/P(n) < Q/4. (26) 

This means that in agreement with Ting and Barnett [3], a cubic medium is completely 
auxetic if Si 2 > and Q > 0, and nonauxetic if 5 , 12 < 0, Q < 0. If none of these pairs of 
inequalities hold, v(m, n) is negative for some pairs (m, n) and positive for other. 



V. CHARACTERISTICS OF ELASTICITY FOR TRANSVERSELY ISOTROPIC 
MEDIA 

Denote the direction of the symmetry axis by c and a component of it by q (i = 1, 2, 3). 
Using the vector c, one can construct two basic second rank tensors, namely a dyad p = c®c 
with components (p) 4 . = pij = CiCj and q = (I 2 — p). Tensors p and q are idempotent, i.e. 
P 2 = P> °L = Q- Their products vanish, p • q = q • p = 0. 

In the case of transversely isotropic media, the tensorial basis consists of five tensors Ei, 
E 2 , E s , F, and G H 

( E i)i^ = E w = VaVku (E 2 ) iifcl = E\f kl = qij q kl /2, 

Gijkl = Gijkl = (PikQjl + PilQjk + PjkQil + PjlQik) /2, 

= (jPijQki + qijPki) I (27) 

7 



3/ ijkl 



+ (E 4 ) 



ijkl 





Ei 


E2 


E3 


E 4 


Ei 


Ei 





E3 





E2 





E2 





E 4 


E3 





E3 





Ei 


E 4 


E 4 





E2 






TABLE I: Multiplication table for E's tensors 

Fijki = Fijki = (qikQji + QjkQu — QijQki) /2. (28) 

Tensors Ej obey the multiplication table HJ 

Similarly as Ei and E2, tensors F and G are idempotent, i.e. F 2 = F, G 2 = G and 

F G = G F = 0, F • E a = E a F = 0, G • E a = E a G = (a = 1, . . . 4). (29) 

For transversely isotropic media (symmetry classes 00, 00/m, oo2, 00m and 00 /mm and 6, 
6, 6/m, 622, Gmm, 6m2 and 6/mmm), the tensor of elastic constants C can be expanded 
in the introduced basis 



C = C1E1 + c 2 E 2 + c s E s + cpF + c G G, 



where 22 



Ci — PijCijkiPki — C33, c 2 — -QijCijkiqki — (Cn + C12) ; 

c s = {PijCijkiqki + QijCijkiPki) = V2C13, 
cf = -FijkiCijki = (Cn — C12) , cq = -GijkiCijki = 2C 44 . 
In the same way we find the coefficients of expansion of S 

si = 5*33, s 2 = (5*11 + S12), s s = \/2S 13 , s F = (Sn - S 12 ) , s G = S M /2. 
Since S is inverse of C and I 4 = (Ei + E 2 + F + G), we find 

S = siEi + s 2 E 2 + s s E s + s F F + s G G = s±, 

where 



c 2 



S\ — , S2 — , S s — , Sp — C F , S G — C G , 



(30) 



se 



se 



se 



(31) 



(32) 



(33) 



(34) 



where s F = [c\Ci — c 2 s ). These equations are in agreement with relation obtained by Boas 
and Mackenzie (cf. Nye 



With the help of the expansion formulas 

1 



22 



2 \f*2i 
J = —Ei H — E? H Ec, 

333 

K = -Ei + -E 2 - — E s + F + G, 

3 3 3 



(35) 



and Eq. (fTTj) . one can check that when C44 = (Cn — C12) /2, C33 = Cn, C13 = C12, 
transversely isotropic media are equivalent to isotropic media. 

Having the explicit form of expansion of S tensor with the help of definitions ®-© for 
transversely isotropic media we obtain 



E-\n) = si (nc) 4 + - (s 2 + s F ) 



1 - (nc) 

2 



^(m,n) 2 2 . 1 / 

■— — — = (si + s F - 2s G ) (mc) (nc) + - (s 2 - s F 
rL{n) 2 



[V2s s + 2s G ) (nc 
1 - 



1 



-s s {(mc) 



1 - (nc) 



1 — (nc 
1 

+ (nc) 2 1 — (mc) 2 I 



mc 



nc 



(36) 



1 



4G(mn) 



S\ + 7rS 2 — V2s s — sq) (mc) 2 (nc) 2 + ]-s F 1 — (mc) 2 1 — (nc)' 



+- 9 sa{( 



mc 



1 - (nc) 



+ nc 



1 — (mc) 2 I 



VI. ELASTIC CHARACTERISTICS OF TETRAGONAL ELASTIC MEDIA 

In the case of tetragonal elastic media, one deals with two groups of classes. To a more 
symmetrical belong 4mm, 42m, 422 and 4/mmm. The remaining less symmetric classes are 
4, 4 and 4/m. 

As the crystallographic directions of all classes of the tetragonal system are mutually 
perpendicular, the related three unit vectors a, b, c also are mutually perpendicular, hence 

a<g>a + b<g>b + c<g)c = I 2 . (37) 

Vector c defines principal axis of symmetry. 
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A. Classes 4mm, 42m, 422 and 4/mmm 



If vectors m and n are perpendicular, 

(ma)(na) + (mb)(nb) + (mc)(nc) = 0. (38) 
Components of tensors E 1; E 2 , E s and G can be expressed by a, b and c j^] 

E ijh = CiCjC k ci, E^ kl = - (a^ + kbj) (a k a t + b k b{) , 

E ijll = ~^p£ iC 3 (°* a l + bkb ^ ' E ifkl = ~/| (°i a J + & A) CkC h E s = ( E 3 + E 4 ) , 

G ijk i = - [{cidj + me,) {c k ai + a k d) + (c^,- + hcj) {c k bi + b k d)} . 

Introduce two further tensors F x and F 2 with components 

F ijki = \ ( a i b j + h i a j) ( a kh + b k ai) , F^ ] kl = ^ (a^ - fyfy) (a k ai - b k b t ) . (39) 

Using the relation of completness (J37j) . one may show that (Fi + F 2 ) = F, where F is defined 
byEq. fify. 

The introduced tensors make up the decomposition of I 4 

I 4 = Ex + E 2 + Fi + F 2 + G. (40) 

Tensors E Q (a = 1, ... ,4) obey multiplication table |U F 1; F 2 , and G are each idempotent, 
and FiF 2 = F 2 Fi = 0, F^G = GF^ = (/? = 1, 2). Products of tensors F^ and E Q vanish. 



As previously, both tensors C and S can represented in form of linear combinations 



22 



C = c^Ed) + 4 2) E 2 + 4 S) E S + 4 ) F 1 + 4 2) F 2 + c G G, (41) 
S = 4 1] E! + 4 2) E 2 + 4 S) E S + s F ~ ) F l + 4 2) F 2 + s G G } 

where Cp (previously we omitted the index E) and c G are the same as in the case of 
transversely isotropic media , whereas Sp (a = 1,2, s) and s G obey Eqs. (|54~jl . The 
coefficients Cp\ Cp and Sp\ Sp are equal 

cg } = 2C 66 , cf = (Cu - C 12 ) , s£ } = S m /2 = 4 2) = (S n - S 12 ) = l/cf. (42) 

If 5*66/2 = (Cu — Cu), then Sp = Sp = sf and Sp 'Fx + SpF 2 = sfF, hence s t i = s±. 
Analogous relations hold for C. 
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Since the first of equations holds, and 



L = Fi + G, M = -Ei + -E 2 - — E. 



(43) 



3 3 3 

when S33 = Six, Sx3 = S12, and 5*66 = 5*44 the compliance tensors of cubic and tetragonal 
higher symmetry media coincide. 

The definitions of Poissson's ratio, Young's and shear moduli (|T|)-(|3*|), and the relation 
(J3HJ) i n the case of tetragonal media of higher symmetry (labeled by the index tl), lead to 
the expressions 



z/fi(m,n) 
E tl (n) 



v m, n 



E(n) 



+ 2s (ma) (mb) (na) (nb) 



(2) 



+ ~ 2 s F 



:ma) 2 - (mb)' 



na 



E^(n) = E~\n) 



+ 24 1) (na) 2 (nb) 2 + ^g ) 



na 



(nb) s 



(44) 



4G t i(m, n) 



1 



4G(m, n) 



+ -sfi [(ma) (nb) + (mb) (na)] 2 

1 2 

+ -s F [(ma) (na) — (mb) (nb)] , 



where 



v m, n 



E(n) 



[sx - 2sq) (mc) 2 (nc) 2 + ^s 2 1 - (mc) 2 



nc 



1 



V2" 



mc 



1 - (nc) 



+ (nc) 2 1 - (mc) 2 I , 



E-Hn) 



,4 1 



sx (nc) + -s 2 



1 - (nc) 



4G(m, n 

4*{( 



a/2s s + 2s G ) (nc 



'i + l;S2 - V2i 



1 - (nc) 



, s G I (mc) 2 (nc) 2 



mc 



nc 



+ (nc) 2 1 — (mc) 2 I . 



B. Classes 4, 4, 422 and 4/m 



In the case of tetragonal less symmetric media (we shall label them by t2) to tensors E a 
[a = 1, . . . , 4), Fi, F 2 , and G, one should add two tensors F 3 and F 4 

F 3 = ^(a<g>b + b<g> a)<g>(a ®a-b<g>b),F 4 = ^(a<g>a-b<g> b)<g>(a <8> b + b <g> a) . (45) 
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Together with tensors Fi, F2 (|39|). they form a subalgebra which is isomorphic to that of 
E a (a = 1, . . . , 4) by having the multiplication table H] with the kernel letter E replaced by 
F. The tensor is obtained when either of F 3 and F 4 multiplies (from the left or the right) 
G or any one of E. 

As previously, the tensors C and S can be constructed as the linear combination 



C = Cl E! + c 2 E 2 + c s E s + c i p ) F 1 + c^F 2 + c^F 8 + c G G, 
S = S «E! + 4 2) E 2 + 4 S) EL S) + 4 1} Fi + 4 2) F 2 + s { ; ] F s + s G G, 



,(2)l 



»1 



(46) 



where F, 



(F 3 + F 4 ), and 



C F 



~w _ J 4 ) _ or* e 1 -*' 

Cp — Cp — 6 F 



„(3) _ „(4) 



5*16) 



and (a = 1, 2, s) are given respectively by Eqs. (J3T|) and (|32p 
In a similar manner to that explained previously, the scalar expressions 



s (i) 

b F±2 



.(2) 



^F ^F 



(2) 
2 ' S F,i2 



(1) 



^F ^F 



(s) 
2 ' S F,t2 



C 



Jl) r (2) 
^F ^F 



2 ■ 



(47) 



(4* 



can be obtained. Using the explixit expressions for coefficients of expansions ()47j1 . we get 

C*16 



5, 



66 



, IP 11 — >->12j — j *->16 

Sp Sp 



Sp 

( Q ) „„J „(") 



(49) 



where Sp = [(Cu — C12) Cqq — 2C 2 6 ]. Expressions relating c ( p' and ; , (a = 1, 2, s) as well 
as Co, and sq, are the same as for tl media. 

Young's, shear, and the Poisson's coefficients for tetragonal media of lower symmetry are 
equal 



v(m, n) 




z/(m, n) 


[ E(n) 


t2 





J tl 



+ 4 s) {(ma)(mb) 
+ (na)(nb) 



^(n) 



E- 1 ^)] +28$ (na) (nb) 



ma 



na 



na) 2 — (nb) 2 
' — b) 2 I 



; mi 



[nbf 



(50) 



1 

4GWm, n) 



4G t i(m,n) 



+ Sp [(ma) (nb) + (mb) (na)] 



J ti 



x [(ma) (na) — (mb) (nb)] 



where 



u(m,n) 



F(n) 



ti 



» i E l ( n ))'tn and 



1 



4G(m,n) 
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can be obtained from 



i/ t i(m,n) 



and 4Gtl ( m n ) as a result of the changes of coefficients Sp to s ( pj 2 (a = 1, 2). 



» 



F t i(n) 
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VII. CONCLUSIONS 



We demonstrated the effectivness of application of tensorial bases of [[V 2 ] 2 ] for derivation 
the explicit expressions for Young's and shear moduli and for Poisson's ratio. These formulae 
are well suited for studying the anisotropy properties of Y, G, and v j^. We considered 
materials of high and middle crystalline symmetry. Work on deriviation of the suitable 
expressions for low symmetry materials and 2D systems is currently in progress. 
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